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Abstract 

The Rayleigh- Taylor instability is a key process in many fields of 
Physics ranging from astrophysics to inertial confinement fusion. It 
is usually analyzed deriving the linearized fluid equations, but the 
physics behind the instability is not always clear. Recent works on 
this instability allow for an very intuitive understanding of the phe- 
nomenon and for a straightforward calculation of the linear growth 
rate. In this Letter, it is shown that the same reasoning allows for 
a direct derivation of the relativistic expression of the linear growth 
rate for an incompressible fluid. 
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1 Introduction 



"What I cannot create, I do not understand" was once found written on 
Richard Feynman's blackboard. Albert Einstein stated that "You do not 
really understand something unless you can explain it to your grandmother." 
The point made by these two great minds was that understanding something 
in physics means you come to the point when it seems obvious and you no 
longer need the equations to derive the result. Among the physical ideas that 
have been examined over the years, the stability /instability concept has been 
extremely fruitful, although its understanding in a given setting would not 
always win Feynman's or Einstein's approval. The stability of a ball inside 
a bowl is very intuitive, and its oscillation frequency when removed from its 
equilibrium position can be calculated just looking at a sketch of the system. 
The instability of a pencil on its tip is equally obvious and frequently cited 
when introducing the concept of an unstable system. In plasma physics. 
Fried could analyze the filamentation instability of two counter-propagating 
particle beams from the very understanding of the physical mechanism at 
work (Fried 1959). Yet, a similar intuitive derivation for the well-known 
two-stream instability is still lacking. 

The Rayleigh- Taylor instability (RTI) plays a key role in many fields 
of physics, and its behavior in connection with inertial confinement fusion 
(ICF) (Pomraning 1990, Kawata et al. 1993), z-pinch physics (Douglas et al. 
2001) or metallic hydrogen generation experiments (Piriz, Lopez Cela, Ser- 
ena Moreno, Tahir & Hoffmann 2006, Lopez Cela et al. 2006) has been the 
topic of many recent works. The RTI occurs when a heavy fluids is accel- 
erated against a light fluid (see Figure [1]). In ICF, a spherical Deuterium- 
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Tritium target is compressed by a Laser. The laser ablates the target, cre- 
ating a low density ablating plasma outside the pellet. During the early 
phase of the compression, the interface between the compressed target (the 
heavy fluid) and the low-density ablating plasma (the light fluid) accelerates 
(Piriz et al. 1997). An observer "sitting" on the interface would then feel a 
force pushing him from the heavy fluid to the light one, resulting in the RT 
unstable configuration pictured on the Figure. In astrophysics, the RTI is 
frequently invoked to explain the filamentary structure of the Crab nebula 
for example (Hester et al. 1996). As the supernova remnant (the dense fiuid) 
decelerates through the interstellar medium (the light fiuid), the interface 
between both is again RT unstable as it experiences an acceleration from the 
heavy to the light medium. 

A Feynman/Einstein like heuristic approach to the RTI (Rayleigh 1900, 
Taylor 1950) was recently provided by Piriz et. al. (Piriz, Cortazar, Lopez 
Cela & Tahir 2006) for an incompressible fiuid. In the usual, normal modes 
approach, where the fiuid equations for both fiuids are linearized, the linear 
growth rate is derived but the basic mechanisms at works remain hidden be- 
hind the equations (Chandrasekhar 1961). In contrast, Piriz et. al. proposed 
a direct derivation of the linear growth rate from the very description of the 
physics involved, "short-cutting" much of the equations. 

2 Non-relativistic approach 

Suppose the interface represented on Figure [T] is initially in equilibrium, both 
incompressible fiuids exerting a pressure upon it. The system is acceler- 
ated downward with an intensity g m/s^. The interface is then displaced by 
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Figure 1: The Rayleigh- Taylor instability. At equilibrium, the pressure on 
both side of the interface is Pq. The pressure variation when moving the 
interface shows the perturbation is unstable if pi > p2- 

z along a distance ^ where k mimics here the wave-number introduced 
in the normal modes approach. The pressure of the upper-fluid at the inter- 
face increases by an amount pigz, while the pressure of the lower fluid also 
increases, but by a quantity p2gz. The upper fluid now pushes the interface 
downward with the pressure po + PiQZ-, and the lower fluid pushes upward 
with the pressure po + P2gz. It is obvious that if po + PiQZ > Vo + P2gz, 
i.e. pi > p2, the pressure balance amplifies the perturbation. Note that 
according to a similar analysis, moving the interface upward equally triggers 
an instability. 

The calculation of the linear growth rate is straightforward from this 
stage. Let us consider a transverse direction, say y, to Fig. [1] so as to 
account for the 3D nature of the system. The surface of the interface over a 
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depth D along the y axis \s S ^ D /k. The force acting upon it thus reads, 

F {pi- P2)gzS = {pi- p2)gz^. (1) 

As it moves, the interface also moves a layer of fluid on both sides. The 
volume of fluid displaced over the surface S should be proportional to S 
itself, and the height of the layer moved on both side is proportional to 
1/k (which is why the normal mode calculation assumes the fluid thickness 
is much larger than 1/k). We can thus write the expression of the mass 
involved in the displacement, 

M oc piS/k + p2S/k = ^^4r^D. (2) 

k-^ 

Since this amount of fluid is displaced by a distance we can write Newton's 
law F = Ma from Eqs.([I]l2]) as. 



where. 



D pi + p2 ^<Pz 2 d'^z 



^2 Pi - P2 J ... 

6 = ■ gk, 4 

Pi + P2 

and a proportionality coefficient equal to unity has been assumed for the 
mass in Eq. ([2]). Equation has exponentially growing solutions if 5^ > 0, 
and 6 is exactly the linear growth rate of the RTI (Rayleigh 1900) where 
the Atwood number At = (pi — P2)/(pi + P2) is immediately identified. It 
is clear now that the sum of the densities relates to the total amount of 
fluid involved in the motion, while the difference relates to the pressure shift 
generated by the displacement. Note that a mass factor different from unity 
in Eq. ([2]) would yield a slightly different denominator for the linear growth 
rate. Indeed, such intuitive calculations frequently yield the correct scalings 
with some pre- factor close to unity. In the present case, the result is exact. 



3 Relativistic version 



The relativistic version of the RTI is especially relevant to Supernovae and 
Gamma- Ray-Bursts physics (Waxman & Piran 1994, Levinsona 2010), where 
ultra-relativistic inhomogeneous flows are involved. Adapting the normal 
modes method to such settings, Allen & Hugues (Allen & Hughes 1984) 
found the relativistic counterpart of Eq. (jl]) , 

Let us now analyze the problem from the intuitive standpoint explained 
above. Starting from Eq. (|4]), where exactly shall we have to introduce rela- 
tivistic expressions? The displacement itself is not relativistic. The interface 
corrugation is a still, initial condition. The only relativistic modification will 
have to do with the inertia of the fluid displaced. A volume of fluid dV has 
the mass pdV in the non- relativistic limit. If the particles it is made of have 
relativistic motion, for example a temperature T such as ksT ~ mc^, the en- 
ergy pdV adds up to the mass within an amount cx pdV/c^. Relativistic fluid 
theory shows indeed that the correct factor is 3 so that the relativistic mass 
density is p + 3p/c^ (Landau & Lifschitz 19876). The density term in Eq. 
([2]) needs therefore to account for this extra inertia. While the interface has 
not been displaced, the pressure is the same on both side and the correction 
per unit of volume reads 3po/c^ for both fluids. Updating the pressure here 
for the corrugated interface would introduce a second order term in z, which 
is neglected in the present linear regime. The relativistic counterpart of Eq. 
([2]) is thus readily obtained replacing pi + p2 by pi + 3po/c^ + P2 + 3pq/c^, 
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and the new linear growth rate is, 



-Sk. (6) 



Qpo/c'^ + Pi + P2 

The calculation starting from the linearized relativistic fluid equations (Allen 
& Hughes 1984) yields the same expression with a term Spo/c^ instead of 



4 Discussion 

Equation ([6]) thus give, up to a numerical factor, the correct value of the 
RTI linear growth rate. We now discuss the discrepancy between the factors 
6 and 8. To do so, we can start from the relativistic Euler equation in the 
absence of gravitational fleld (Landau & Lifschitz 1987a), 

(^ + e)u'^ = ^-uu'^ (7) 
^P^^)"" dxi^ dx^ dx>^' ^ ' 

where p is the pressure, e the energy density, = (ct, r) and = (7, 7v/c). 
As previously said, our problem is relativistic in the sense that the energy 
density can be so, not by virtue of some relativistic velocity of the fluid 
elements. Setting thus 7 = 1, one can check that the temporal component 
{i = 0) of the equation above yields v ■ Vp = 0. The spatial part (i = 1, 2, 3) 
then gives 

Y dp „ p + e f d \ , , 

Neglecting v ■ Vv as a second order quantity and adding the acceleration 
gives the premise of Allen & Hugues' Eq. (3), 



Equations fl7f^ show that that the correct relativistic inertia is not simply 
the energy density e, but the energy density plus the pressure, e +p. Setting 
then e = pc^ + 3p gives p + e = pc^ + Ap for both fluids, from which the factor 
8 eventually arises. 
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